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IMPORTANT NOTICE 


Texas Instruments (Tl) reserves the right to make changes in the 
devices or the device specifications identified in this publication 
without notice. Tl advises its customers to obtain the latest version 
of device specifications to verify, before placing orders, that the 
information being relied upon by the customer is current. 


Tl warrants performance of its semiconductor products, including SNJ 
and SMJ devices, to current specifications in accordance with Tl’s 
standard warranty. Testing and other quality control techniques are 
utilized to the extent Tl deems such testing necessary to support this 
warranty. Unless mandated by government requirements, specific 
testing of all parameters of each device is not necessarily performed. 


In the absence of written agreement to the contrary, Tl assumes no 
liability for Tl applications assistance, customer’s product design, or 
infringement of patents or copyrights of third parties by or arising from 
use of semiconductor devices described herein. Nor does TI warrant 
or represent that any license, either express or implied, is granted 
under any patent right, copyright, or other intellectual property right 
of Tl covering or relating to any combination, machine, or process in 
which such semiconductor device might be or are used. 


Copyright © 1985, Texas Instruments Incorporated 


INTRODUCTION 


The TMS32020 Digital Signal Processor is a fixed-point 
16/32-bit microprocessor. However, it can also perform 
floating-point computations at a speed comparable to some 
dedicated floating-point processors. 

The purpose of this application report is to analyze an 
implementation of floating-point addition, multiplication, and 
division on the TMS32020. The floating-point single- 
precision standard proposed by the IEEE will be examined. 
Using this standard, the TMS32020 performs a floating-point 
multiplication in 7.8 microseconds, a floating-point addition 
in 15.4 microseconds, and a floating-point division in 22.8 
microseconds. 

To illustrate floating-point formats and the tradeoffs 
involved in making a choice between different floating-point 
formats, a review of floating-point arithmetic notation and 
of addition, multiplication, and division algorithms is first 
presented. 


FLOATING-POINT NOTATION 


The floating-point number f may be written in floating- 
point format as 


f = mxbe 
where 
m = mantissa 
b = base 
e = exponent 


For example, 6,789,320 may be written as 
0.6789320 x 10/7 


In this case, 


m 0.6789320 
b = 10 
e = 7 


The two floating-point numbers f; and f7 may be written as 


fj m, xbel 
fp = my xbe2 


Floating-point addition/subtraction, multiplication, and 
division for fj and f2 are defined as follows: 


fj Xf) = (my+my xb ©1—©2)) xbel if ey =e0 (1) 
Or 
= (mj x b~ €2-eD) +m) x be2 if e, <e? 
f) xf) = my xm xb@1 + &2) (2) 


(m4/m3) xb (1~€2) (3) 


f/f 


A cursory examination of these expressions reveals 
some of the factors involved in the implementation of 
floating-point arithmetic. For addition, it is necessary to shift 
the mantissa of the floating-point number which has the 
smaller exponent to the right by the difference in the 
magnitude of the two exponents. This is shown in the 
multiplication by the terms 


b~ (€1—-€2) and b~ (€2—e1) 


This right shift can result in mantissa underflow. There 
are also possibilities for mantissa overflow. Addition and 
subtraction of exponents can lead to exponent underflow and 
overflow. To alleviate underflow and overflow, it is 
necessary to decide on some scheme for roundoff. For a 
detailed description and analysis of underflow and overflow 
conditions and rounding schemes, see reference 1. 

It is desirable to have all numbers normalized, i.e., the 
mantissas of f; and fz have the most significant digit in the 
leftmost position. This provides the representation with the 
greatest accuracy possible for a fixed mantissa length. The 
result of any floating-point operation must also be 
normalized. The factors associated with normalization, 
overflow, and other characteristics of floating-point 
implementations are best illustrated with a few examples. 

Consider the addition of two binary floating-point 
numbers f; and f7 where 


0.10100 x 2011 
0.11100 x 2001 


fy 
f2 


Both of these numbers are normalized, i.e., the first 
bit after the binary point is a 1. Addition requires equal 
exponents, so the fractions are aligned by shifting right the 
one with the smaller exponent and adjusting the smaller 
exponent. This yields 


fo = 0.00111 x2911 
Then, 


fy +f) = 0.10100 x 2011 +0.00111 x 2011 
= 0.11011 x2011 =f, 


The sum may overflow the left end by one digit, thus 
requiring a postaddition adjustment or renormalization step. 
Since it is assumed that the register is only of a finite length, 
this renormalization will result in the loss of the lowest order 
bit. 

Another example illustrates the overflow past the most 
significant bit. With an assumed register length of five, let 


0.11100 x 2011 
0.10101 x 2001 


fy 
f2 


Then, 


0.11100 x20ll=f, 
+ 0.0010101 x2011 =f, 
1.0000101 x 2011 =f, 


The significance of the two digits underlined in the right 
part of the mantissa is suspect, since it is assumed that the 
corresponding bits of fj are zero. The left underlined digit 
is the overflow past the most significant bit. To finish the 
addition, f3 is shifted to the right and the exponent adjusted 
accordingly. Thus, 


1.0000101 x 2011 =f, 


The shift of the fraction and the adjustment of the exponent 
yield 


0.10000101 x 2100 = fz 
The result may be rounded, giving 
0.10001 x 2100 = fz 
or truncated, giving 
0.10000 x 2100 = f, 
FLOATING-POINT ALGORITHMS 


Multiplication Algorithm 

The algorithm for normalized floating-point 
multiplication is illustrated in Figure 1. This algorithm is an 
implementation of Equation 2 in the section on floating-point 
notation. The floating-point numbers being multiplied are A 
and B written as 


A = ma Xb*A and B = mp xb‘B 
The result is 
C = mcxbsc 
For the resulting mc, there are three special cases. The 


mc may be zero, in which case there is a branch to Step 
10 to set C=O. If mc #0, then the most significant bit will 


MA mB 


4 MULTIPLY MANTISSAS 
mc =m, X mB 


CA eB 
ADD EXPONENTS 
2 ec = ea + 6B 
3 TEST FOR SPECIAL CASES OF mc 


NORMALIZED 





TEST FOR SPECIAL CASES OF ec 


OVERFLOW 
OR 
UNDERFLOW 


10 | SET SPECIAL VALUES OF RESULT 






C=AxB 


Figure 1. Floating-Point Multiplication 


be in either the first or second leftmost bit. If the most 
significant bit is in the second leftmost bit, then a left shift 
of mc is necessary (see Step 5). Otherwise, C is already in 
normalized form, and there is a branch to Step 6. 

In Step 6, the desired rounding scheme is implemented. 
After this rounding, it is possible that mc will overflow (see 
Step 7). In this case, it is necessary to right-shift mc one 
bit (see Step 8). Special cases of ec , are tested for in Step 9. 


If there is an overflow or underflow of ec, it is corrected 
in Step 10. Otherwise, the result is in range, and the 
calculation is complete. 


Addition Algorithm 


The implementation of normalized floating-point. 


addition is more involved than for multiplication. This 
addition algorithm, outlined in Figure 2, is an implementation 
of Equation 1 in the section on floating-point notation. 
In Step 1, eq and ep are compared to determine ec. 
For this illustration of the algorithm, it is assumed that 
eA sep. The right shift (d) required to align ma, is 
determined in Step 2. The procedure in Step 3 implements 
the right shift of ma. In Step 4, the extra bits of ma are 
discarded by using the desired rounding technique. The 
mantissas of A and B are then added in Step 5. 


A B 


| [emcee Tac RO Sone 
: amc eons oe 
s [Risrmanrens oar na are 
a 
s [eowarasne ne oa 
+ [earns core 


NOR 
OVERFLOW eee 









k LEADING ZEROS 


RIGHT SHIFT mc ONE BIT LEFT SHIFT mc k BITS 
ec =ec + 1 ec = ec - k 
FORCE EXPONENT 
9 ec = 0 










DISPOSE OF EXTRA BITS: 
ROUNDING OR TRUNCATION 
TEST FOR OVERFLOW OF mc 


NO OVERFLOW OVERFLOW 


RIGHT SHIFT mc ONE BIT 
ec =ec + 1 
TEST FOR SPECIAL CASES OF ec 


IN RANGE 


12 





OVERFLOW 
OR 
UNDERFLOW 


14 | SET SPECIAL VALUES OF RESULT 


C=A+8B 





Figure 2. Floating-Point Addition 


Now, the procedure becomes somewhat more involv- 
ed. The mc may be zero, in which case there is a branch to 
Step 9 which sets ec = 0; a branch to Step 14 sets the special 
value of the result. The mc may overflow, making a right 
shift of one necessary (see Step 7). The mc may have k 
leading zeroes; therefore, a left shift of k is required. This 
normalization step is generally the most involved and time- 
consuming step to perform. The procedures in Steps 10, 11, 
and 12 round the mc, test for a possible overflow due to 
the rounding, and adjust ec accordingly. The special case 
of ec is determined in Step 13. Finally, after Step 14, the 
sum C = A + Bis formed. 


Division Agorithm 

Floating-point division is more sophisticated than 
multiplication and addition since fixed-point processors such 
as the TMS32020 are not inherently capable of performing 
division. For example, 1/3 = 0.3333...; only an approx- 
imation can be calculated since 1/3 must be represented in 
a finite number of terms. Several algorithms can be im- 
plemented to find good approximations of such numbers. The 
algorithm implemented in this report is shown in Figure 3. 

Step 1 shows the equivalent of A/B. In Step 2, the latter 
term is expanded using a power series of 1/(1 + X), where 
€ (BLO/BHI) is X (€ simply denotes that the term is right- 
shifted 16 bits forming the least significant bits of a 32-bit 
number). The third term in the power series only affects the 
LSB of a 32-bit result; therefore, this term and all the 
following terms can be dropped, as shown in Step 3. 

The equation in Step 3 can be implemented on the 
TMS32020 in two steps. Assuming that the result is a 32-bit 
number Q and that it is composed of a 16-bit QHI and a 16-bit 
QLO, think of the equation in Step 3 in the following 
manner: A/B = Q — eX. The first term is a fair approx- 
imation of the result Q, and the second term is a correction 
term to obtain a better approximation. With this in mind, 
it can be shown that (AHI + eALO)/BHI will give a 16-bit 
quotient and a 16-bit remainder. Due to the architecture of 
the TMS32020, the 16-bit quotient will be in the low word 
of the accumulator and the remainder will be in the high word 
of the accumulator after the division. Since it is desirable 


A divided by B 


AHI + €ALO 
BH! + €BLO 
1 ; 1 
QWORDSIZE 316 


where A 


wo 
tol oll 








AHI + €ALO = AHI +: € 
STEP: “BHeELO” — 
1+ € (a5 ae [BOY ) 
BHI 
AHI + E€ALO BLO BLO \2 
STEP 2: a! eee Te + €2 i 
BHI BHI BHI 


AHI + AHI + €ALO BLO AHI + €ALO 
STEP 3: ~ 2 e( — ) ( = ) 


BHI Hl 


Figure 3. Division Equation 


to have a floating-point result, the remainder must be divid- 
ed by BHI to obtain the low word of the quotient. Now QHI 
and QLO have been calculated. When placing Q into the cor- 
rection term (equation in Step 3), note that Q is equal to QHI 
+ QLO. It can be shown that QLO will have no effect on 
the result since the correction term is multiplied by e. 
Therefore, to calculate A divided by B, simply implement 
the following equation: 
A A BLO 


— = —— —e(—— x qui} 
B CBHI BHI 


where the division is fixed binary (left-shifts and subtracts). 

Figure 4 shows the implementation of the division 
algorithm that was outlined in Figure 3. 

In Step 1, the dividend is right-shifted four times to 
prevent an overflow. Note that the result is not shifted left 
to compensate for this shift, because the normalization routine 
automatically does this. The shift causes the dividend to be 
limited to 27 significant bits instead of 31. In Step 2, a binary 
divide (left-shifts and subtracts) is implemented on the 
dividend by the high 16 bits of the divisor. The 32-bit result 
contains a quotient in the low 16 bits of the accumulator, 
and a remainder (R1) in the high 16 bits of the accumulator. 
R1 is left-shifted fifteen places in Step 3. The new RI is 
divided by BHI in Step 4 to calculate the lower 16 bits of 
the quotient. 

The quotient has now been approximated. The 32-bit 
result is composed of QHI and QLO, as shown in Figure 3. 
To obtain a better approximation, one term in the power 
series expansion must be added to the quotient. Therefore, 
the procedure in Step 5 calculates a 16-bit correction term, 
which is then added (or subtracted since it is the term 
following the “‘1’’ in the power series) to the 32-bit quotient. 

Testing for an overflow of the resulting mantissa is 
necessary. Since the dividend was left-shifted four places, 
the resulting quotient will not be negative if an overflow 
occurred. To detect an overflow, bit 28 in the quotient must 
be tested. If this bit is a 1, an overflow occurred; if it is a 
Q, no overflow occurred. If an overflow has occurred, the 
exponent must be incremented. Finally, it is necessary to 
normalize the quotient and output the results. 


A DIVIDED BY B 


7 


_ 
oO 
_= 
oT 
oO 


WHERE A= ALO 


31 
SHIFT “‘A’’ RIGHT FOUR TIMES TO 
1 PROTECT FROM OVERFLOW. 


A/BHI = 32-BIT RESULT. 
2 HIGH 16 BITS ARE REMAINDER #1 (R1). 


LOW 16 BITS ARE HIGH QUOTIENT (QHI). 





SHIFT R1 LEFT FIFTEEN TIMES. 
EQUIVALENT TO R1 x 219. 


R1/BHI = 32-BIT RESULT. 
4 HIGH 16 BITS ARE REMAINDER #2 (R2). 
LOW 16 BITS ARE LOW QUOTIENT (QLO). 


5 MULTIPLY QHI BY BLO AND DIVIDE BY BHI. 
(QHI x BLO)/BHI = CORRECTION TERM. 


SUBTRACT 16-BIT CORRECTION TERM 
6 FROM 32-BIT QUOTIENT. 
(QHI | QLO) — (0 | CT) = RESULT. 





NO OVERFLOW 


Cexp = Aexp ~ Bexp 


OVERFLOW 


Cexp = Aexp - Bexp + 1 
9 NORMALIZE RESULT. 


Figure 4. Floating-Point Division 


IEEE FLOATING-POINT 
SINGLE-PRECISION FORMAT 


Of interest is a set of formats known as the IEEE 
standard. This IEEE recommended format consists of a 
variety of precision formats (single, double, single-extended, 
and double-extended). The IEEE has also proposed several 
techniques for handling special cases such as overflow, 
underflow, + oo, and rounding. For complete details, the 
reader is referred to the proposed IEEE standard.2 

The single-precision format is a 32-bit format consisting 
of a 1-bit sign field s, an 8-bit biased exponent e, and a 23-bit 
fraction f (see Figure 5). The value of a binary floating-point 
number X is determined as follows: 


X = (-1)8 x 2(¢7-127) x Lf 


31 30 23 22 0 


Figure 5. IEEE Floating-Point Single-Precision Format 


The advantage of this format is that it is structured in 
such a way as to provide easy storage and straightforward 
input/output operations on 8-, 16- and 32-bit processors. The 
disadvantage with this format is that the large mantissa will 
generally span several words of memory. 


FLOATING-POINT IMPLEMENTATION 


IEEE Implementation 

The IEEE single-precision format is described here as 
it applies to the addition, multiplication, and division 
algorithms. In these floating-point routines written for the 
TMS32020, all results are truncated to 31 bits to provide 
more flexibility in the user’s development of a rounding 
scheme suitable for his application. The representations of 
+ oo are ignored so that the user can decide how to handle 
these exceptions in a manner that is appropriate for his 
particular application. 


I/O Considerations 

The first consideration is the internal representation of 
the binary floating-point number. If the number is read into 
the TMS32020 as two 16-bit words, some processing is then 
necessary to put the floating-point number into a 
representation which is easier to process. The representation 
used in the TMS32020 programs in the appendices is shown 
in Figure 6. This internal representation may be arrived at 
by a simple manipulation of the IEEE bit fields. For this 
particular algorithm, it is assumed that the floating-point 
number is input to the TMS32020 as the four 16-bit fields 
shown in Figure 6. However, the user can easily supply his 


own routine to arrive at this format from two 16-bit inputs 
to the TMS32020 where the inputs contain the IEEE single- 
precision format. 

The format in Figure 6 was chosen to minimize the 
execution time of the floating-point addition, multiplication, 
and division routines. The format of the result is shown in 
Figure 7. Notice that it is identical to the format in Figure 
5 except for CLO. CLO has its 16 most significant bits valid 
for both the addition, multiplication, and division routines. 


Normalization 

Since the floating-point routines require normalization, 
a partial binary search algorithm is implemented in the 
addition and division routines in the appendices. To begin 
the normalization routine, note that all mantissas can be 
considered to be positive with the format used for the result 
shown in Figure 7. The binary search for the most significant 
bit (the leftmost 1 since the mantissa is positive) is illustrated 
in Figure 8. 

The first move is to split the result into CHI and CLO. 
If CHI + 0, the most significant bit (MSB) is the CHI; 
otherwise, it is the CLO. For this example, it is in CLO. 


15 141312 1110 9 8 7 6 5 4 3 2 1 90 
ASIGN 
(0 IF POSITIVE, 
OR — 1 IF NEGATIVE) 


15 14131211109 8 7 6 5 43 2 1 #90 


AHI 
(NORMALIZED) 


0 $1 f (most significant 14 bits) 
15 1413 12 1110 9 8 7 6 5 4 3 2 1 #90 


ALO Of (least significant 9 bits) 0 0 0 0 0 90 


15 1413121110 9 8 7 6 5 4 3 2 1 #9 


AEXP 


Figure 6. Floating-Point Representation 


15 141312 1110 9 8 7 6 5 4 3 2 1 «0 


CSIGN 

(0 OR —1) 
15 14 13 12 1110 9 8 7 6 5 4 3 2 1 =«0 

CHI 

(NORMALIZED) 


15 14 13 12 1110 9 8 7 6 5 4 3 2 1 #0 


CLO 


15 14131211109 8 7 6 5 4 3 2 1 #0 


CEXP 


Figure 7. Result Representation 


31 


C 


0 


0000000 0000000000001 010031100071 71 +1 


31 CHI 16 


000000000000000 0 


31 16 


0001010011 000%1 1 ~=41 


BLZ NOFLOW 
RPT TTEEN 
NORM 


0101003110001 1100 


y 


15 CLO | 0 


0001010011000 71 1 ~=1 


15 0 


0000000000000 0 0 0 


GOTO NOFLOW ON OVERFLOW. 
TTEEN = 13, PERFORM 14 ‘‘NORM”’. 


1 Vis 0 


0000000000000 00 0 


Figure 8. Partial Binary Search 


The next step is to form a 32-bit result with CLO in the most 
significant word position. It is now possible for the MSB 
to be in the highest bit location since CLO has been left- 
shifted 16 times. If this is the case, an overflow has occurred, 
and the result must be right-shifted once. The normalization 
routine tests this by branching to NOFLOW if the result is 
negative. If the number is not negative, the normalization 
can continue. 

The NORM instruction is used in the repeat mode to 
complete the normalization. Note that this whole 
normalization routine can be replaced by the following two 
instructions: RPTK 29 and NORM. The RPTK instruction 
causes the NORM instruction to be repeated 30 times, thus 
normalizing a 32-bit number. This method is not 
implemented here due to the timing. These two instructions 
always take 31 cycles to normalize a 32-bit number. The 
normalization routine here takes only 22 cycles (worst case) 
for normalizing a 32-bit number. Therefore, if program space 
is More important than timing efficiency, it is best to replace 
the normalization routine with these two instructions. 


Added Precision 
As illustrated in Figure 7, the 16 most significant bits 
of CLO are valid, i.e., C is valid for 31 places beyond the 


6 


binary point. Oftentimes the user is not as concerned with 
the IEEE standard as in being certain that he has enough 
accuracy for his particular application. Since the TMS32020 
uses 16-bit words, the routines in the appendices implicitly 
maintain a 30-bit mantissa. They also implicitly use a 16-bit 
exponent. If the user desires this added accuracy and dynamic 
range, then it is readily implementable with no additional 
cost in execution time. The normalization for the addition, 


' aS mentioned previously, operates over the entire 32-bit 


accumulator. For the strict IEEE format, the user will only 
want to normalize over the 25 most significant bits of the 
accumulator. The structure of the normalization routine 
makes this modification simple. 

The routines in the appendices make no provision for 
the representation of + oo and exponent underflow and 
overflow. The user of the routines should consider the degree 
of significance of these results and the way they should be 
handled for his particular application. Since these routines 
are written to operate at maximum speed, truncation of results 
is used. If the user desires to implement a rounding scheme, 
then he will also need to check for the possibility of overflow 
due to the rounding scheme. This step is shown in the 
multiplication, addition, and division flowcharts (see Figures 
1, 2, and 3). 


SUMMARY 


The TMS32020 may be used to perform floating-point 
operations with great accuracy, wide dynamic range, and 
high-speed execution. The design engineer has the 
responsibility of deciding what type of floating-point format 
is best for his application. To aid in understanding floating- 
point operations, several examples have been given that 
illustrate the manipulations necessary to implement floating- 
point addition, multiplication, and division algorithms. 
Flowcharts for these algorithms are also included. The 
appendices contain the TMS32020 code for the IEEE 
floating-point single-precision format used in addition, 
multiplication, and division. The addition and multiplication 
routines may also be used without modification to implement 
a format with up to a 30-bit mantissa and a 16-bit exponent 
without any increase in execution time. 
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Qi, Lor Bl AHI (NORMALT ZED) 


tere cnane ovene wneee eupey autor GEES cunge euED ees Sane: SEEPS emote etees Sere SOEED soees sonen 


i i 


eeaes eames eeee creed Soca: stem HBuee epEDe suns cout amwee count epeee artes mpeee cnnes muame mene 


ry FRITS 1-0; ALD 


eenem qneme canes conse GanEh TERRE artes cates menat eopty cOnED out cnegh eoEED ene conse anees teas 


cowed ewwee ennne scene secs seers etewe qtaey cneme cum sent enue seems Cneee woere mpmeE eneee wevet 


AEAP (=127- 70 2283 


cores aenme CUeed epeee GmnOe eOEES Omeee wHESE cEEED aDttD Somes mntES tocms cncen Aeewe erwee enaee moses 


OOS 
OOF? 
OOSE 
OO 
OO30 
OOS 
OOs2 
OOs3 
OOn4 
ORE 


NOAA 


TO CORRESPOND WITH IEEE FORMAT, 
PIPL Cig otf “ae Seo se apd) 
INSTEAD OF 1.F # 2 ##E, AND SUBTRACT Lev’ FROM E. 


THE FINAL FORMAT TS THE SAME AS THE INITIAL FORMAT 
EXCEPT THAT FOR CLO WE HAVE: 


ceeme meres ewene wosee comme ooeet arene cnmam cones aunts CURSE SCRE! EeEED cores comms SunED coner muEer 


1G. VL le CLO 


meee 


OO37 ALL 14 BITS OF CLO ARE VALID. ANYTHING FAST THESE HAS 
OOS BREEN TRUNCATED, 


me ee 


Oo4o PEHPRRHALKKAPHHLHHHHAAERREHSHERAEREH EEE EERE ERE EER ERE EHHEEY 
OoO4] ¥ : 
OO4s 
OoO4s 
Oodd 
OOAs 
O44, 
OO4? PEUPUPAKPHHKHHEH ERE HHSH AREAS ALESHAHHRSEHHERSHHAEHHHHEERS 
OAs # 

OO47 OOO Acres 
OOO O00 ASTON Ertl 

OOe 1 OOO, AEKF Ela! 

OOS OOS AHI Er 


WORST CASE CEXCLUDING INITIALIZATION AND T/a): 
eeac MICROSECONDS, 

THIS TIMING INCLUDES THE NORMALIZATION, 

WORDS OF PROGRAM MEMORY: 2 


me OK 


ae ody 
ee a ee ee 


x 


a a 
“Senet 


mm TD pf ib po etd 


we? 


Aes OOO8 ALO Pill 
OoOe4 OOO4 BSTON EC 


QO QO0O5 BEX EGIL 
Oo Oa048 HHI Fill 


C-2 


NOSTDT 


OOE7 
Ooss 
Oo? 
OO40 
OO4G1 
ONG 
O04 
OO4A4 
OO6gS 
OOS 
OO4? 
OO4S 
OO4 
OO7O 
OO7 1 
OO72 
OO7 3 
0074 
007% 
OO74 
OO? 7 
OO7S 
OO7F 
OO8O 
OOF 1 


OOS. 


ee 


mm ght, tm. ete 


OMS 
OOS 
OOS 


"sf 
mth 


OC 
OO87 
OOSS 
OOD 
OVO 


OCD | 
Oo? 


2 
o 


cn 
OOE 
OODVA 
OOP? 
OOS 
OCI 
OVO 
Oo 
OLG2 
OLos 
44d. 
O1oOs 
OL Oe, 
O17 
OLcns 


OLay 


OOOO 
OOO1 
OOF 


OOO 
‘eleler-! 
OOOE 
OOD 
OO 7 


OOOE 
OO0F 
OO0A 
OOOB 
OOO 


QOOOTI 
OOOE 
OOOF 
QO LO 
QOL] 
Oo LT? 
OOS 
OO1 4 
OOS 
COOLS 
OGL? 
OOS 
OOo. 
OOLA 
MOOR 
OO 
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OGO7 
OOGE 
OOO? 
OOO 
OOOR 
OOO: 
MOOT 


QOOOE 


OOOF 
O010 
Oo. 
OOS 
OOS 
OO 14 
OOS 
OOA 
O017 
OOS 
OO? 


eer 
CEO? 


SS 


Died 
OF00 
CBO? 
SOAQ 


ae 


$ 
ce 
cweat Be 


~ 
a] 
ae 


Lad 
naa 
FFFF 
AOOK 
Doo, 


19000 
4014 
CAQ4 
AQQE 
AGT 
4015 
PAGS 
ANIL 
CAGE 
6017 
CALE 
4018 
LAO 
é01F 
ADO 
AOC? 


BLO 

MS 1GN 
mEXF 
HI 
mame 
NEWGONE 
TEMF 
Four 
iM 

iL 

Fil 

Re 

CL 
MiO00 
CINE 
THREE 
FITEEN 
THIRTY 
TTEEN 


a a a 


Frill 
Fru 
eri 
Er 
Frit 
Frat 
Ertl 
Ertl 
Er 
Eri 
Eri 
Eri 
ELL 
Emu 
Et 
Eri 
Emu 
Emil 
Eril 


1 ow. 


J 
‘ore 


Leneobe owe 2% 
+ i ot obi 
foneete oo 


12 
13 
14 


f 


« 
ove! 


1é 
A 
if 
1? 


20 


— = 


whee. sual 


INITIALIZATION 


LOPE 
=aXM 
LARF 
LALE 


RE TE 
IN 

LARF 
LARE. 
LAL 


SACL 
LALE 


SACL 
AE 
SACL 
LACE 
SACL 
LACE 
SACL 
LACK 
SAL 
LACE 


4 


1 
ARL, 2z00 


7 
H+ FAG 
a 
ARO, 6G 
>FFFF 


NEGONE 
21000 


MLO 
4 

Fife 

1 

CINE 
THREE 
cs 
PLITEEN 
“4 
THIRTY 
ti 
TTEEN 


EXP 


PCO, 7 4. 345 


FAGE OO002 


BEGIN ON FAGE 4, 
SET &IGN EXTENSION, 


READ NUMBERS INTO BLOCK Bo, 


MLEAR EXPONENT REGISTER. 


NEGONE = 1 


MiOooo = =1ooo 


FOUR = 4 


MINE = 1 


THREE = 3 


Pa PRES 1S 


THIRTY = 20 


TTEEN = 12 


LEAR CEXF 
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Olid 
OLd1 
OL1e 
OLis 
O114 
OL15 
O114 


O17 
Ois 
O11? 
OL2o 
O1e1 
Ole 
Ole3 
O124 
Olf3 
O1246 
Ole? 
O25 
Oe? 
0130 
Osi 
G1se 
O1S35 
O14 
Oe 
OVS4 
Q4 37 
Oss 
GUY 
O140 
O14) 
0142 
O45 
O144 
O145 
O1 44 
0147 


O1 48 
o14y 
O50 
O1S1 
OVS2 
OVS 
O154 
OLS 
O1LE4 
O37 
OYSs 


ee 


0140 
O41 
O142 


OO 


OOLF 
OOLO 
OO 1 


OO2s 


OOS 
OSA 
OOF 


OOZL 


One 7 
OoOFs 
OO? 
002A 
OOZE 
OO 
OO2T 
OSE 
OOF 


N30 
OO 
O08? 
OO 


OOs4 


OOSS 
OOS 4 
OOS7 


OOS 
OOR? 
OOS 
OOSR 
OOS 
OosD 
OOSE 
OOS 
OOad 
O04 1 
0047 
OO4S 
0044 


O04% 
OO44, 
O47 
QOO4S 
OOdy 
OO4A 


F000 
1oo04 
FASC 
Oo2s 


2OOL 


4no7 
AIS 
4AB14 
Cel? 


4B17 
ATO: 
6214 
AQOF 
SF ii 
4E17 
A704 
4212 


4010 


OF 
a? 
CE14 


» 4817 


470O& 
&é(Q) 4% 
4Q00F 
4210 
1O13 
6008 
S508 
200A 
4E1i4 
FeeO 
OO" f 
LOLS 
AOD 
BOO] 
1005 
Ooo? 
400% 


200A 
FARO 
OO4E 
4008 
4°08 
4Bi+ 


xs 


a oe 
.’ 
% 
ae 


NOV 


me 


NORMAL 


FOGE ooo 


FINISHED WITH INITIALIZATION 


LAC 
SACL 
SLB 
BZ 


LAC 
SACL 
ZALH 
ADDS 
RET 
SFR 


RET 
IBC: 
S ACH 
SACL 
LAL 
RET 
SR: 
SAH 
AL 


pal 
MFY 
FAL 
RET 
SUBIC 
SACL. 
ZALH 
ADT 
SUB 
Sac. 
SACH 
LAC 
ANT 
BR? 


LAC 
SACL 
LAC 
SLB 
ADD 
SACL 


ASTON 
mS TGN 
BS TGN 
tik 


NEGONE 
ms TGN 
AHI 
ALU 
THREE 


FITEEN 
BHI 

Rt 

cM 
Riis 
FITEEN 
BHI 

Re 

CL 


LM 
BLO 


FITEEN 
BHI 
ick 

iM 

GL 

ae 
Cho 
CHI 
CHI 
MiaQo00 
NOOVE 


CINE 
EAP 
AEXP 
RE XF 
ILE XF 
EXP 


NORMALIZE 


LAL 
BZ 


ZALH 
AnH 
RPT 


CHI 
LU 


HI 
ML 
TIEEN 


CSTGN = ASSIGN, IF ASTON = BSTGN, 


ELSE, CSIGN = -1. 


SHIFT DIVIDEND TO FROTECT FROM OVERFLOW. 


MM = AHI; ALO / BHI, Ri = REMAINDER. 
HIGH ACCUMULATOR RETAINS REMAINDER. 


f BA 
(Ri # Se#175) 


(Ri # 2##15) 
COMPUTES 


WIVES GL, AND Re. 
/ BHI. 


HIGH ACCUMULATOR RETAINS REMAINDER. 


MORRECTION TERM = (QM # BLO) 
MOMPUTES (iM # BLO). 


/ BHI. 


COMPUTES (iM # BLO) / BHI. 


BMG — OFCL = CHIICLO 


DIO AN OVERFLOW OCCURS 


TF NOT, GOTO NOOVF. 
ELSE, INCREMENT CEXP. 


IOMPUTE RESULTING EXPONENT. 


DOES CHI HAVE THE MSB? 


Ip YES, NURMALIZE RESULT. 


WILL PERFORM 14 "NORMS", 


C-3_ 


NOS TOT 


OLAS 


O144 


O1E5 
OQ144 


O1A47 
O16 
O169 


OL70 
O171 
O172 
O73 
0174 
O175 
OL74 
O177 
OL7#8 
O17? 
O10 
O1S1 
0152 
O1ES 
o134 
o1es 
OLS4 
O187 
o1ss 


o1e9 


OO4B 
OO4C 
OoO4T 
QOO4E 
OO4F 
OO8O 
OOS 
OOS 
OOS 


OOF 4 


OQEE 


OO 


O07 
OOS 
oOas? 
OoO5A 
OOSB 


OOM 


NO ERRORS , 


32020 FAMILY MACRO ASSEMBLER PITQ, 7 84, 348 TaaS5317 O3-27-85 
##e# PRERELEASE ### 


MEAS 
FREQ 
OOR7 
400B Ll 
Faso 
HORS 
4Bi¥ 
TEAS 
FREQ 
0057 


me oe oe ok ok Ok 


ITEO& NOFLOW 
CEI 


i i A a 


+ 
OUITRLET 


Cer 
NO WARNINGS 


PAGE o0n4 


NORM 
EB CWITRPLUT GO GLITFUT RESULTS, 


ZALH PLO HERE, IF CLO HAS MSE. 
Bee NOFLOW DIQ BIT SEARCH CALISE OVERFLOW? 


RFT TTEEN IF NOT, NORMALIZE RESULT. 
NORM 
B CUT RUT GO OUTPUT RESULT. 


FINISHED WITH NORMALIZATION 
OVERFLOW OCCURRED DURING BIT SEARCH 


RSXM ROAM FOR LOICAL RIGHT SHIFT. 
=FR PERFORM RIGA? Siri 


TARE CARE OF EXPONENT & NORMALIZED MANTIS3A, 
THEN OUTPUT RESULTS. 


SACH THT SAVE NORMALIZED MANTISSA€. 
SACL CLO 

LARF 1 RESET POINTER, 

RPT THREE MUTPUT RESULTS, CSIGN, CEXP, CHI, AND CLO. 
CLIT #+, FAQ 

LOGE WAIT FOR INTERRUPT. 
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